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We investigate the properties of the pion quasiparticle in the low-temperature phase of 
two-flavor QCD on the lattice with support from chiral effective theory. We find that the pion 
quasiparticle mass is significantly reduced compared to its value in the vacuum, by contrast 
with the static screening mass, which increases with temperature. By a simple argument, 
near the chiral limit the two masses are expected to determine the quasiparticle dispersion 
relation. Analyzing two-point functions of the axial charge density at non-vanishing spatial 
momentum, we find that the predicted dispersion relation and the residue of the pion pole 
are simultaneously consistent with the lattice data at low momentum. The test, based on 
fits to the correlation functions, is conhrmed by a second analysis using the Backus-Gilbert 
method. 


I. INTRODUCTION 


Identifying the spectrum of excitations of strongly interacting matter at finite temperature is 
of central importance to understanding the nature of the medium. These excitations are encoded 
as poles in thermal correlation functions. For some quantum numbers, the existence of weakly 
interacting probes coupling to local operators make it possible to measure the properties of these 
excitations experimentally, at least in principle. The prime example of such a probe is the photon. 
In practice, a medium is created in heavy-ion collisions which appears to reach thermal equilibrium 
locally, the temperature decreasing with time. Therefore a weighted average of thermal photon or 
dilepton spectra over the spacetime history of the ‘fireball’ is obtained (see e.g. mm and Refs, 
therein). 

In the low-temperature phase, it is natural to ask how close the properties of the excitations 
are to those of the known hadrons at zero temperature. Viewed globally, the spectrum does not 
appear to change much until temperatures close to the transition temperature are reached, where 
the rapid crossover to a deconfined and chirally symmetric phase takes place. This conclusion is 
based on the success of the hadron resonance gas model in describing equilibrium properties of the 
medium (particularly the equation of state and quark number susceptibilities) computed in lattice 
QCD [SHS], and on its success in describing particle yields in heavy-ion collisions [HIT]. However, 
reliable information about individual excitations is sparse. 

Here we extend our study [8] of the pion at finite temperature in two-flavor lattice QCD with 
support from a thermal chiral effective theory [MI]. In El we performed a temperature-scan on 
16 X 32^ ensembles; here we focus on one temperature (T = 170 MeV) on a fine lattice (24 x 64^) 
with high statistics, for which we also have a reference zero-temperature ensemble at the same bare 
parameters. Our study shows that the zero-temperature pion mass ‘splits’ at finite temperature 
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into a lower pion quasiparticle mass and a higher pion screening mass, 

r = 0 : pion mass = 267(2) MeV 

^ \ 

T = 169MeV : quasiparticle mass = 223(4)MeV screening mass = 303(4)MeV. 

While the quasiparticle mass is the real-part of a pole of the retarded correlator |p| = 0) 

of the pseudoscalar density in the frequency variable, the (static) screening mass is a pole of 
Gr{oj = 0, p) in the spatial momentum |p| and represents an inverse spatial correlation length. 
The pion quasiparticle mass can be extracted model-independently near the chiral limit due to 
the dominance of its contribution to the two-point function of the axial charge density and of the 
pseudoscalar density. By contrast with the mass, we find that the decay constant associated with 
the pion quasiparticle practically retains its zero-temperature value. 

By a simple argument, the dispersion relation of the pion at low momenta is given by a single 
parameter u (see Eq. 0) below), which in the chiral limit corresponds to the group velocity of 
the excitation. As in [8] , we determine this parameter as the ratio of the quasiparticle mass to the 
screening mass, u ~ 0.74. As a new aspect, we then test whether the so-determined parameter 
u correctly predicts the momentum-dependence of the pion energy by looking at the two-point 
function of the axial charge density at non-vanishing spatial momentum. An important observation 
is that the chiral Ward identities also predict the residue of the pion pole in the axial-current two- 
point functions. Due to the difficulty of extracting real-time information from Euclidean correlation 
functions, testing simultaneously the predictions for the pole and the residue proves to be essential 
to improve the discriminative power of the analysis. Since the chiral predictions are only expected 
to be valid at sufficiently small momenta, we also provide an estimate of the range of validity of 
the effective theory. 

Our analysis method of lattice correlation functions is based on fits, where the ansatz is mo¬ 
tivated by the chiral effective theory at small frequencies and on perturbation theory at high fre¬ 
quency. We also present an alternative analysis, which starts by generating model-independently 
a locally averaged spectral function by following the Backus-Gilbert inversion method [IMZl. In 
a second step, a pion pole contribution with the predicted dispersion relation is assumed, allowing 
us to obtain an estimate for the residue. The advantage of this alternative analysis is that we 
do not have to formulate an explicit ansatz for the spectral density of the non-pion contributions. 
This point is particularly relevant since at finite spatial momentum, axial-vector excitations do 
contribute to the two-point function of the axial charge density. 

The paper is organized as follows. Section [II] contains an overview of the theory expectations 
concerning the two-point functions of the axial current at finite temperature. Sections |III| and |IV| 
present the lattice QCD calculation, and our conclusions are given in section |Vj In appendix 
the tensor structure of the axial-current two-point functions at finite temperature is given; in 
appendix we derive the contribution of the pion to the four independent tensor structures, 
thus determining all the relevant residues. Einally, appendix contains a table with the lattice 
correlator data. 


II. THEORY BACKGROUND 

We work in the Euclidean path integral formalism, and our notation and conventions follow 
those used in [8]. The vector and axial-vector vectors and the pseudoscalar density are given by 

= ^(x)7/,y = V!(x)7^75y V’(a:), i^“(x) = (1) 
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where is the isospin-doublet quark field. In appendix [Xj we provide a decomposition in momen¬ 
tum space of the Lorentz structure of the two-point functions of the axial current. For a general 
momentum p, they are entirely described by four ‘form factors’, which in the rest frame of the 
thermal medium are functions of po and p^. At zero-temperature, the four functions reduce to 
two functions of p^, one longitudinal and one transverse. The partially-conserved axial current 
(PCAC) relation relates the two-point function (P“(x)P^(0)) of the pseudoscalar density, as well 
as the (A“(x)P^(0)) correlation functions to the aforementioned form factors. 

In this work, we investigate the following static screening correlators, 

T) = j dxo <fx^ (a^(x)A|(0)) , (2) 

6^^ GU^3, T) = I dxo d^xx_ (p“(x)P'’(0)) , (3) 

where x± = (xi,X 2 ). Time-dependent correlators with a general spatial momentum p will also 
play a crucial role. 


5^^Ga{xo,T,p) = ^ 

1 (AS(i).4j(0)) , 

(4) 

h“'Gp(xo,r,p) = ^ 

f cfx (p“(i)P'’(0)) . 

(5) 


They are related by Fourier transformations to the form factors defined in appendix [A] for instance 

G\{x3,T) = j ^e-*^’3-3nL,i(0,p2). (6) 


The correlators Ga{xo,T,0) and G\{x 3 ,T) are only sensitive to the longitudinal form factor 
n^’^; these were the cases considered in [8]. At non-vanishing momentum however, the correlator 
Ga{xo,T,p) is sensitive to three independent form factors II'^ and 11^’^. 

At long distances, the screening correlator G\{x 3 ,T) is given by 

( 7 ) 

which defines the screening pion mass m-,^ and the associated decay constanl|^/.,r- The Gell-Mann- 
Oakes-Renner relation 


flrnl = -miipij;) ( 8 ) 

holds to leading order in the chiral expansion. From Eqs. (§ and the low-momentum analytic 
structure of the longitudinal form factor 11^’^ reads 

nRi(o,p2) = ^^, p^o. (9) 

More generally, expanding the denominator in the frequency, 

+ ml —> p‘^ + ml + + ..., ( 10 ) 


3 The normalization convention is such that at zero temperature ~ 92 MeV. 
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FIG. 1. Trajectory in (w,p)-plane of the pion pole in the pseudoscalar retarded correlator. At negative 
the pole corresponds to the pion screening mass, at positive p^ it corresponds to the pion quasiparticle. The 
slope is the value of u^iT). 


it follows that a quasiparticle (pole in the retarded correlator as a function of frequency) with the 
dispersion relation mm 


UJt 


up = n(T)Vm 2 +p 2 ( 11 ) 

exists at low momentsj^ The remarkable aspect is that the parameter u determines both the (real 
part of the) dispersion relation of the quasiparticle and the ratio of the quasiparticle mass to the 
screening mass. A graphical interpretation of the dual role of the parameter u is given in Fig. 
Here the trajectory in the frequency-momentum plane of a pole in the retarded correlator of the 
pseudoscalar densitj0 corresponds to a static screening state at and to a real-time 

quasiparticle at small positive p^. In [8], the parameter u was determined using lattice correlation 
functions at vanishing spatial momentum via the two estimators 


Um — 


Uf = 


4m2 Gp(xo,T,0) 


ml Ga(xo,T,0) 

fn^TT 

xo=/3/2_ 


2Ga(/ 3/2, T, 0) sinh(u/m7r/3/2) 


( 12 ) 


(13) 


In doing so, the parametric dominance of the pion in the time-dependent Euclidean correlator at 
small quark masses is exploited. Good agreement was found between Uf and Um at T ~ 150 MeV 
for a zero-temperature pion mass of 305 MeV. Any departure of u from unity clearly represents 
a breaking of Lorentz invariance due to thermal effects. In this work, one of our goals is to test 
whether the parameter u determined from the ratio of the quasiparticle to the screening mass, as 
in [18], really does predict the dispersion relation of the quasiparticle, as in Eq. 0- In order to 
carry out this goal, we perform an analysis of the time-dependent Euclidean correlator GAixo, T, p) 
in terms of the spectral function p^. They are related as follows. 


Ga{xo,T,p) 



duo p (w, p) 


cosh(ti;(/3/2 — xq)) 
sinh(io/3/2) 


(14) 


^ In this argument, the imaginary part of the frequency-pole is neglected. A more sophisticated argument is required 
to show that the damping rate of the pion quasiparticle is indeed parametrically subleading US- 
® Recall that the momentum-space Euclidean correlator GB(tu„,p) is related to the retarded correlator via 
= GE(oJn,p) for > 0 [H]. 
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First we recall that at zero temperature, the Lorentz structure of the axial-current two-point 
function implies the following momentum dependence of the pion pole contribution, 


GAixo,T = 0,p) 


|xo|^oo O' ® 




Res(a;p) 


H 


(15) 


where the residue is here given by 

Res(w°) = (/° a;° = (p^ + 

In terms of the spectral functions, this correlator corresponds to 

T = 0, p) = Res(a;°) - ( 0 ;°)^) -h •.. 


(16) 


(17) 


At the other end of the spectrum, in the high frequency region, a leading-order perturbative 
calculation (see for instance [l9]) yields 


p^(a;,r,p) = 6'(w2-4m2-p2)^^(p2-h6m2), w 


00. 


(18) 


We note that at non-zero momentum, the correlator Ga{xo,T = 0,p) receives contributions from 
axial vector mesons via the transverse form factor (see appendix [A]) , 


GA(xo,r = 0,p)= f 


+ P^) + 2^1 + P^) 


LPo + P 


Po + P 


(19) 


The spectral functions associated with the form factors II"'" and 11^ are measured experimentally 
in r decays [20] . Since II"'" describes by itself the two-point function of spatial components of the 
axial current at vanishing spatial momentum, it cannot contain the pion pole. The latter is entirely 
contained in the form factor 11^. 

At finite temperature, the pion pole appears in all three form factors contributing to Ga{xo, T, p). 
Altogether, the pion contribution to the spectral function is predicted to have the form 
p^{uj,T,p) = Res(wp) 5 ( 0 ;^ — LOp), with the dispersion relation given by Eq. (11) and the residue 
by (see appendix [B|) 


Res(a;p) = f^{ml + p^). 

For later use we also define the pion quasiparticle decay constant /* via 

Res(wo) = 


The contribution to the Euclidean correlator then reads 


Ga{xo,T, p) 


Res(a;p) cosh(a;p(/3/2 — xq)) ^ 

2ca;p sinh(a;p/3/2) 


( 20 ) 


( 21 ) 


( 22 ) 


Whether the residue determined through fits to lattice correlation functions agrees with Eq. (20) 
provides a cross-check that the low-energy effective description is working. 
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III. LATTICE SETUP 

In this section we describe the analysis performed on a finite temperature ensemble of size 
24 X 64^ with two degenerate dynamical light flavors. The short direction is interpreted as time 
and therefore the temperature is T = l/24a = 169(3)MeV while the spatial extent amounts to 
L = 64a = 3.1 fm. The fields admit thermal boundary conditions in time and periodic boundary 
conditions in space. We use the Wilson plaquette action [2T] and the 0(a) improved Wilson 
fermion action with a non-perturbatively determined Csw coefficient |22]. The configurations were 
generated using the MP-HMC algorithm [23] following the implementation described in [23] based 
on Liisher’s DD-HMC package [25| . In addition, we use a 128 x 64^, effectively zero temperature 
ensemble that was made available to us through the CLS effort (labelled as 07 in |26|) with all 
bare parameters identical to our finite temperature ensemble. The pion mass takes a value of 
ttItt = 270MeV |26| such that itIt^L = 4.2. This additional zero temperature test ensemble allows 
us to compare thermal observables in a straightforward manner with their corresponding “effective 
zero-temperature” value calculated in the 07 ensemble. 


A. PCAC mass 


In order to check that our thermal ensemble indeed yields the same physical quark mass as its 
corresponding zero-temperature counterpart (07), we use the definition of the quark mass based 
on the PCAC (partially conserved axial current) relation 


l/dxodV(9fP^™P(x)P(0) 

m^CAc{x2.) - - J. (P(x}F(0)} 


where in the improvement process 


X_L = (xi,X2) 


(23) 




(24) 


The derivative 5™’’ is the improved lattice discretized version of the derivative following |29) . The 


non-perturbatively calculated coefficient ca was taken from I30|. Note that since relation ( |23| ) stems 
from an operator identity, we are free to choose the direction of measurement. On the thermal 
ensemble, the spatial directions are longer; therefore, by measuring along these directions we obtain 
a longer plateau and thus more accurate determinations of the PCAC mass. The extraction can be 
carried out by performing a fit to a constant in the range where a plateau is observed. Its central 
value and error, given in Table |T| are in very good agreement with the ones quoted in m- 


B. Pseudoscalar and axial-vector correlators 


Our goal is to calculate the temperature-dependent coefficient u{T) that parametrizes the pion 
dispersion relation 0 . In |S|, we defined two estimators Uf and Um that yielded consistent 
results up to T ~ 170MeV for the case of two degenerate light flavors with = 2GeV) ~ 15 

MeV; at that quark mass, the crossover region is located aronnd Tc ~ 211 MeV [32|. In the 
thermal ensemble we are analyzing here, we have = 2GeV) = 12.8(l)MeV (see Table |^, 

and therefore expect a slightly lower value of the transition temperature. Nevertheless, this should 
not affect the applicability of the chiral expansion around (T, m = 0) with T < Tq, as discussed in 

m 
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x^ja 


FIG. 2. The PCAC masses. The renormalization factors and Zp are included, as well as the conversion 
factor from the SF to MS at a scale of ^ = 2GeV, which amounts to 0.968(20) [2S]. We also plot the result 
along the XQ-direction to show that indeed both are compatible. This can be interpreted as a check that 
cutoff effects are indeed small for this value of the lattice spacing. 



5.50 

K 

0.13671 

Csw 

1.751496 

TN^= 2 i [MeV] 

169(3) 

Tv^=i 28 [MeV] 

32(1) 

a [fm] [2^ 

0.0486(4)(5) 

za Eg 

0.793(4) 

Zp [26] 

0.5184(53) 

fn^^/T{n = 2GeV) 

0.0757(7) 


TABLE I. Summary of the main parameters for the 24 x 64^ finite temperature ensemble as well as for the 
128 X 64^ zero temperature ensemble labelled as 07 in |26j . The quark mass is computed at and normalized 
with the T = l/24a temperature. The statistics collected for two-point functions is respectively 360 and 
149 configurations at N^. = 24 and Nj. = 128, with respectively 64 and 16 point sources per configuration, 
exploiting the translational invariance of the system. 


We use the correlators defined in Eqs. ([2j|^ with the spatial momenta given by 

P = Pn = (0,0, 27rn/L). (25) 

The improvement of the axial current was already introduced in Eq. ( |24[ ). Note that since all 
two-point functions belong to the adjoint (or isovector) representation of SU(A^/) {Nf = 2), the 
contributions of quark disconnected diagrams cancel out. The renormalization program is carried 
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FIG. 3. Effective ‘cosh’ masses for the 07 zero-temperature ensemble in the (Aq^o) channel for n = 0,1. 
The values of w® = are given in units of the temperature T = l/24a corresponding to our thermal 
ensemble. 


out such that 


Ga{xo,T,p) = {ZA{gi)fGA{xo,gQ,T,p) (26) 

Gp{xo, T, p) = {Zp{g^o)fGp{xo, gl T, p) (27) 

and similarly for the screening correlators; the value of the coefficients Za and Zp can be found in 
Table m 


IV. ANALYSIS OF LATTICE TWO-POINT FUNCTIONS 

After the preliminary work presented in the previous section, we turn to the analysis of corre¬ 
lation functions in order to extract the pion properties. 


A. The zero-temperature case 


As a benchmark we analyze zero-temperature data on the 07 ensemble. Here we are able to 
obtain the pion energy Wp by fitting a constant to the effective mass. The pion energy corresponding 
to p = 0 and p = (0,0, 27r/L) can be read off from the plot in Fig. The dominance of the pion 
contribution, particularly in the zero-momentum case, is clearly very strong. Performing a linear fit 
to (cUp)^ as a function of p^, we obtain for the slope v?{T ~ 0) = 1.01(6), consistently with Lorentz 

defined by Eq. (16), indeed turns out to be independent of the 


invariance. The decay constant f'^, 
momentum. 


Once the ground state dominates the correlator Ga{xo, T ~ 0, p), one-state cosh fits of the form 
Ai cosh(a;p(T/2 — xq)) with T = 128a are applied and the results are summarized in Table [n| The 
values for cjq and are in very good agreement with the ones quoted in m- 
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n 

Ai/T^ 


X^/d.o.f 

f^/T 

Res(a;p) 

0 

8.4(3) X 10-3 

1.579(12) 

0.05 

0.599(8) 

0.89(3) 

1 

5.3(4) X 10-4 

2.88(3) 

0.4 

0.629(12) 

3.27(15) 


TABLE 11. Properties of the pion at zero temperature. The index n denotes the momentum p„ induced 
and ojp^ corresponds to the energy of the state (in particular, Wq is the pion mass). All errors are purely 
statistical, and all renormalization factors are included. The fit interval begins at Xg/a = 6 for the zero- 
momentum case and at xo/a = 15 for one unit of momentum in view of the effective mass plot of Fig. 
Dimensionful quantities are normalized with T = l/24a. 


B. The screening quantities and at finite temperature 

A detailed description of how the extraction is carried out can be found in [ 8 ]. Here, we highlight 
the basic relations that lead to the extraction of the screening quantities /tt and m-,^ and therefore 
to the values of the estimators uj and Um- 

• The screening mass m,r is calculated by fitting the correlation function Gp{xs,T) with a 

two-state ansatz of the form Aj cosh[mj(L /2 — X 3 )]) with masses m, and amplitudes 

Ai. The value obtained for the ground-state mass is compatible with the value obtained from 
the ‘cosh’ mass which is defined as the positive root of the following equation (see Fig. [^, 

(?|>(x3,r) _ cosh[mcosh(a:3 + a/2)(x3 - L/2)] ^ 

Gp{x 3 + a, T) cosh[mcosh(iC 3 + a/2){x3 + a- L/2)]' 

• We determine the screening pion decay constant from the correlator G'^(x 3 ,T) via Eq. ([^ 
by applying again a two-state ‘cosh’ ansatz. The screening pion mass ttItt also appears in 
G\. We use this fact as a consistency check, but due to the better signal to noise ratio of 
the pseudoscalar channel, we quote the value extracted from Gp as our final result for m.,^. 


C. Thermal time-dependent correlators at zero momentum 


The estimators uj and Um for the pion velocity u are defined in Eqs. ( T2p(3 ). Apart from the 

PC AC mass, /jr and ttItt, they involve the time-dependent correlators Ga{xo, T, 0) and Gp{xo,T, 0). 

The difference between Uf and Um can be explained as follows. The estimator Um is based on the 

dominance of the pion contribution in the correlators Ga and Gp at xq = 13/2; the estimator Uf is 

based on assuming that the residue is given by the screening quantities (as predicted by the chiral 

effective theory), Res(wo) = and the dominance of the pion contribution in Ga only. The 

dominance in Ga is less strong an assumption than the assumption that the pion dominates Gp, 

2 

since their spectral functions are related by pp{u}) = —^^pa{uj) at zero spatial momentum. We 
summarize results for Uf and Um in Table [111} 

The chiral expansion around (T, m = 0) proposed in |llj assumes that one is sufficiently close to 
the chiral limit. In this limit, the screening pion mass m-Tr vanishes and the coefficient u{T) is indeed 
the velocity of a massless pion quasiparticle in the presence of a thermal bath. A deviation from 
unity corresponds to a violation of boost invariance. At finite but small quark mass, we showed 
in [ 8 ] that the consistency of uj and Um serves as an indicator for the applicability of the chiral 
expansion. Based on the results of Table III we conclude that they are indeed consistent. The 


results, in particular for Uf, are in good agreement with the results obtained in [ 8 ] on ensembles 
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FIG. 4. Effective ‘cosh’ mass plot for the screening mass m.^-. The plateau has been chosen to begin at 
xs/a = 19. 



1.79(2) 

UIT 

0.46(1) 

w/ 

0.76(1) 

Um 

0.74(1) 

Uf/Ujn 

1.02(1) 

uJo/T 

1.32(2) 

fJT 

0.62(1) 

Res(a;o)/T4 

0.68(2) 


TABLE III. Summary of the results for the AA = 24 thermal ensemble. All renormalization factors are 
included and the errors are purely statistical. The value of ujq is calculated using ujq = UmixiTr. In the same 
way fl = fTrjum- The value of the residue is obtained according to Eq. (20), Res(a;o) = 


with a coarser lattice spacing and a slightly heavier quark mass. Note that to leading order, u is 
expected to be independent of the quark mass. 

The ‘reconstructed’ correlator is defined as the thermal Euclidean correlator that would 
be realized if the spectral function remained the zero-temperature one. We compute it following 
the method first proposed in [33] . Figure shows the difference between the thermal correlator 
and the reconstructed correlator. There is a statistically significant difference between the two 
correlators, which shows that a change must take place in the spectral function. Because the 
difference is very weakly dependent on time, the change must take place in the low-frequency part 
of the spectral function. We expect from the thermal chiral effective theory that the change is due 
to a modification of the mass and/or the residue of the pion quasiparticle. Using the numbers of 
Table III, the changes amount respectively to 


^ = 0.836(14), 


Wn 


jk = l-03(2). 


(29) 













11 


0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 
0 

FIG. 5. The difference AGyi(|p|,Xo) = [GA(a;o,T, p) — GJ|‘^(a:o, T, p)]/T^ of the thermal correlator and the 
reconstructed correlator for |p| = 0 and 27r/L. 


AGA(0,a;o) 

AGA{2Tr/L,xo) 


10 


11 


12 


xo/a 


We thus observe that while the pion decay constant remains unchanged at the precision level of a 
few percent, the pion mass decreases by 16%. Qualitatively, these results are consistent with the 
two-loop results in zero-temperature chiral perturbation theory given in |34[ [35]. Future lattice 
calculations approaching the chiral limit would allow for a quantitative comparison. 


D. The spectral function at non-zero momentum 


As the next step, we test the functional form of Eq. © at non-zero momentum. The relevant 
real-time pion states with energy Wp have a non-zero overlap with the operator f d?x e*P^j4o(x); 
furthermore, the spectral function pA becomes independent of u) in the ultraviolet, rather than 
growing like We therefore expect to have the best sensitivity to the pion contribution in the 
correlator G/i(xo,T, p). 

At finite temperature, the analysis of the correlator Gyi(xo,T, p) is more involved than at zero 
temperature: only at sufficiently small quark masses and momenta, and not too small xq is the 
correlator parametrically dominated by the pion pole. Therefore we proceed by formulating a fit 
ansatz to take into account the non-pion contributions. The combination of Eqs. ( |17| ) and (18) 
motivates an ansatz for the spectral function reading 

/5'^(w,p) = Ai(p)sinh(a;^/2)d(a;-u;p)-hA2(p)^^ (^1 - ^(a; - c). (30) 

The corresponding form of the correlation function then reads 


G^(xo, T, p) = Ai(p) cosh(u;p(/3/2 - xq)) ^ 2 (p) 


Nc 

247r2 


g-cxo g-c(P-xo) 


Xo 


/3 - Xo 


(31) 


We fit Ga{xq,T, p) with the ansatz given in Eq. (31) for the momenta pn = (0,0, 27rn/L) with 
n = 1, 2,3,4, 5. The fit interval is chosen to be xo/a G [5,12] in order to avoid cutoff effects. There 
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are four parameters involved, j4i(p), Wp, ^ 2 (p) and c. Leaving Wp as a fit parameter led to poorly 
constrained fits. Therefore we set the value of ojp to the prediction of Eq. 0 in order to test 
whether the data can be described in this way. Motivated by the expected large-w behavior of 
the spectral function, we quote the rescaled parameter A 2 = ■ Note that the quark mass is 

negligible compared to all the non-vanishing |p| values considered here. The expected value of A^ 


is of order unity, in view of Eq. (18). Eq. (22) allows us to establish the relation between the fit 


parameter Ai(p) and the residue itself, 

Res(a;p) = 274i(p)a;p sinh(u;p/3/2). 

One may further convert the result for the residue into a parameter 6(p), defined by 

Res(wp) = fl{ml + p2)(l + 5(p)). 


(32) 


(33) 


Erom the chiral prediction Eq. 
effective description is working. The results are summarized in Table IV 


we thus expect 6(p) to be small compared to unity if the 


n 

Ax/T^ 

Wp„/T 

A 2 

c/T 

Res(wp^)/T4 

b 

X^/d.o.f 

1 

2.95(4) X 10-1 

2.19(3) 

1.78(8) 

6.7(3) 

1.72(6) 

-0.08(3) 

0.06 

2 

1.40(5) X 10-1 

3.73(6) 

1.26(2) 

6.1(1) 

3.3(2) 

-0.39(4) 

0.15 

3 

4.9(3) X 10-2 

5.40(9) 

1.19(1) 

7.7(1) 

3.9(5) 

-0.65(4) 

0.35 

4 

1.7(2) X 10-2 

7.1(1) 

1.15(1) 

9.67(9) 

4.21(7) 

-0.78(3) 

0.49 

5 

4(1) X 10-3 

8.8(1) 

1.12(1) 

11.7(1) 

3(1) 

-0.89(3) 

1.04 


TABLE IV. Results of fits to the axial-charge density correlator at non-vanishing momentum p„. All errors 
quoted are statistical, and all renormalization factors are included. The quantity Wp/T is not a fit parameter; 
rather it is set to the value predicted by Eq. (111 with u[T) = Um = 0.74(1). 


The fits provide a good description of the data; see the y^/d.o.f values and Eig. We observe 
that at the smallest momentum, |p| ~ 400 MeV, b{p) really is small, pointing to a successful check 
of the chiral prediction. At higher momenta, the negative, order unity value of 6(p) indicates that 
the residue of the pion pole is reduced. It should also be remembered that at higher momenta, 
neglecting the width of the quasiparticle is bound to be an increasingly poor approximation. The 
coefficient A 2 is expected to be of order unity from the treelevel prediction Eq. (18). Indeed the 


numbers in Table |IV| are of order unity. One reason for the relatively large value of the coefficient 
at the smallest momentum could be that axial-vector excitations are contributing around the 
threshold c, thus adding spectral weight. The value of the threshold at |p| ~ 400 MeV, is about 
l.lGeV, a value we consider to be reasonable given that the mass of the lightest axial-vector meson 


in nature is ~ 1.2GeV. 


In order to gauge the discriminative power of the test, it is interesting to ask whether a rather 
different model is consistent with the lattice data on G'yi(xo, p, T). We assume for this purpose 
that the dispersion relation and the residue have the same p-dependence as at zero temperature. 


p = y/ojQ + p^, and obtain for n = 1 an equally good description of the data. 


We therefore set ^ ^ _ 

with a value of the residue Res(a;p) = 3.01(4) not too different from (/*)^(u;q -|-p^) = 2.84(7). The 
other fit parameters take the values A 2 = 2.42(17) and c/T = 10.2(4). While the perturbative 
coefficient and the threshold values seem less plausible to us, we cannot completely exclude this 
model on the basis of the lattice data. 

To summarize, we have found that the dispersion relation of the pion quasiparticle is consistent 
with Eq. (0, the parameter u being determined at vanishing spatial momentum. In order to test 
the dependence of our results on the fit ansatz made, in the next section we apply the Backus- 
Gilbert method. 
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FIG. 6. Correlation functions G^(a:o, T, p„)/T^ with p„ = (0,0, 27rn/L), together with the fits resulting 
from the four-parameter Ansatz of Eq. The corresponding parameter values are given in Table |IV| All 
renormalization constants are included. 


E. The Backus-Gilbert method for p^(w,p) 


The Backus-Gilbert method is a method suitable for inverting integral equations like Eq. (14). It 
has been studied in many contexts (see e.g. [IMZ!). While it has not been applied in lattice QCD, 
to our knowledge, the central notion of resolution function was used in [36]. We first describe the 


method in some generality. It is a completely model-independent approach since no ansatz needs 
to be made for the spectral function. 

The goal is to solve the integral equation 


G{xi) = / dujf{u})K{xi,uj) Xj / 0 Vi 

Jo 


(34) 


for the unknown function /(w), given the kernel K{xi,Lo) and given data on G{xi). The idea is to 
define an estimator /(w) 


POO 

f{ui)= / 6{O,0j)f{uj)du) 

Jo 


(35) 


where <5(a), a;) is called the resolution function or averaging kernel. It is a smooth function concen¬ 
trated around some reference value a), normalized according to duj6{Cd, w) = 1, and parametrized 
at fixed u) by coefficients qi{Cd), 


i 


(36) 


so that / is obtained according to 


f{^) = '^G{xi)qi{u). 
i=l 


( 37 ) 
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The goal is then to minimize the width of the resolution function. Minimizing the second moment 
of 6{u),uj) in its second argument around its first argument yields 


EjWr\u)R{x,) 

Ek,im)w^iHo)R{xiy 


(38) 


where 



du}K{xi, 


R{xi) 



K{xi,uj)duj. 


ujyK{xj,ui), 


(39) 

(40) 


We remark that /(ca) equals /(w) if the latter is constant. 

The matrix Wij{u) is very close to being singular. This is the reason why, when trying to use 
a data set with error bars, one needs to regulate and change the matrix W to 






XWij + (1 - A)5, 




0 < A < 1 


(41) 


where Sij is the covariance matrix of the data. The value of A controls the trade-off between 
resolution and stability. For values of A close to 1, we obtain the best possible resolution. However 
the results tend to be unstable since the matrix is poorly conditioned and large cancellations take 
place among the terms in Eq. (37). Reducing A improves the stability of the result at the cost of 
deteriorating the frequency resolution. 



FIG. 7. Left: Some examples of resolution functions for different values of A centered at w/T. Right: 
Estimators /3(a;,p„)/T^ for n = 0,1,2,3 together with its error shown as a band. The vertical colored 
dashed lines correspond to the locations of the expected positions of the poles according to Eq. (11) 
with u{T) = Um- The black horizontal lines correspond to the treelevel asymptotic values of All 

renormalization constants have been taken into account as well as the improvement program on the axial 
correlators. Dimensionful quantities have been made dimensionless by the appropriate power of T = l/24o. 


We apply this method to Eq. (14). In order to regularize the finite-temperature kernel at w = 0 
we rewrite the equation as 


Ga(xo,T,p) 



f ^ f cosh(a;(/3/2 - xo)) A 

\tanh(u;/2 )) \ cosh(a;/3/2) ) 

' -V-" 

=K(xo,u)) 


0 


(42) 
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This defines our estimator p at uj = co, 


p{u},p) = 


f 


' 1 ^ /^ Pa{oj,p) 

i t,nh(a,?/2) 


(43) 


After regulating the problem via the covariance matrix Sij as in Eq. (41), the inversion is carried out 


via Singular Value Decomposition. This offers the opportunity to diagnose how badly conditioned 
the matrix is. With all quantities made dimensionless by applying appropriate powers of the 
temperature, we choose A = 2 • 10“^ in the following. Typical condition numbers of the regularized 
matrix in Eq. 0 are ~ 10®. The situation gets worse when A approaches unity, as explained 
above. The results for zero momentum and the first three units of momentum are shown in the 
right panel of Fig. As in the case of the fit, we included the points of the correlation in the 
interval xo/a = [5,12] so Wij{0) is a n x n symmetric matrix with n = 8. The value of A was 
chosen equal to 0.002. This yields a relative error on p of ~ 3 — 5%, while the resolution function is 
displayed in the left panel of Fig.[^ One direct observation is the fact that the expected asymptotic 
behavior for large values of w is reproduced very well. 



FIG. 8. The effective residue Res(wp,a;)BG as defined in Eq. (441. Left: No momentum induced, p = 0. 
Right: One unit of momentum induced, pi = (0,0, 27r/L). The grey band is the expectation in terms of 
screening quantities, Eq. ( [20| . All renormalization factors are included. The errors arise from the statistical 
uncertainty. The values of Wp are indicated by dashed vertical lines. 


The right panel of Fig. [^also shows the expected positions of the poles that follow from Eq. ( |11[ ) 
as vertical colored dashed lines. We now want to test the p dependence of the residue Res(a;p) via 
the following argument. If we assume that, for a given value of w, p{u}, p) is dominated by the pion 
pole, we obtain the following estimator for the residue. 


Res(wp,a;)BG = 


2 wp tanh(wp/3/2)/>(a;, p) 
(5(a;, Wp) 


(44) 


Here we treat iOp as input and calculate it using Eq. 0 with the value of u = Um = 0.74(1) 
determined at zero momentum. The result as a function of co is shown in Fig. [^for zero and one 
unit of momentum. The natural choice where Res(a;p, a;)BG is expected to be the best estimator 
of the residue is at a; ss cOp. Looking at Fig. Isl one sees that approximately around this value the 


curve intercepts the grey band, which represents the prediction Eq. (20). The latter is particularly 
well verified at zero momentum, while the agreement at jpj = jpij ~ 400MeV is at the ten percent 
level. These observations provide a further test that the pion dispersion relation Eq. © predicted 




















16 


by the thermal chiral effective theory is consistent with the lattice two-point function of the axial 
charge density. 

Comparing the method followed in this subsection with the previous method based on a global fit 
to the spectral function, the former has the advantage of not requiring an explicit parametrization 
of the non-pion contributions to the spectral function. This observation may be useful in other 
lattice studies of spectral functions. 


V. CONCLUSION 

We have found that the pion quasiparticle mass is reduced significantly by thermal effects 
compared to its vacuum value - unlike the pion screening mass, which increases. Also, the energy 
cost of giving the pion quasiparticle a momentum is significantly reduced, since the ‘velocity’ is 
well below unity, u k, 0.74. We have tested that the pion indeed admits a modified dispersion 
relation, Eq. by analyzing lattice two-point functions. The test is based on requiring the 

consistency with the lattice data of the combined chiral prediction for the dispersion relation and 
the residue of the pion pole in the two-point function of the axial-charge density. These conclusions 
could be strengthened further by repeating the calculation at smaller quark masses and with higher 
statistics. Discretization errors should also be studied. Having a higher resolution in momentum 
could help in assessing the region of validity of the chiral effective theory. 

Assuming the results hold to further scrutiny, one may wonder how much a modified mass and 
dispersion relation of the pion affects (a) the freeze-out mechanism in heavy-ion collisions and (b) 
the predictions of the hadron resonance gas model for equilibrium properties. In answering the 
latter question, one must take into account that the modification of the pion dispersion relation 
is due to the presence of hadrons in the medium, and issues of double counting arise. However, 
the following estimates may provide a useful first idea of the size of the effect. At the temperature 
of 169MeV in the two-flavor theory that we have been discussing, with a zero-temperature pion 
mass of 270MeV, we estimate, using the hadron resonance gas model, an isovector quark number 
susceptibilit}|^ amounting to XsjT'^ = 0.42. In the HRG model, the pion contribute^ Xs/7"^|pions = 
4/3^ f ( 27 rp /b(p)(1 + /b(p)) = 0.28. If the spatial-momentum integral in the pion contribution is 
cut off at pmax = 400 MeV (roughly the range of validity of the chiral effective theory that we 
found), the contribution is reduced to 0.11. If we instead use the modified dispersion relation with 
the lower quasiparticle mass ujo = 223MeV and u = 0.74, the contribution for p < Pmax amounts 
again to 0.28. It is unclear whether one should include a contribution from higher momenta, given 
that the thermal width of the pion may then not be negligible. The numbers above illustrate that 
the contribution of the pion to the quark number susceptibility might not be as strongly affected 
as one may at first think. However, the contribution comes from softer pions, which implies a 
reduced amplitude of the transport peak in the two-point function of the vector current V)“(x), an 
effect that can be tested in lattice simulations. 

Determining the dispersion relation of a non-Goldstone hadron would be interesting to see 
whether the relatively strong change we have seen in the pion properties is specific to chiral dy¬ 
namics. In general, a kinetic theory description allows one to use as primary degrees of freedom 
the quasiparticles specific to the temperature of interest. It is therefore much broader in applica¬ 
bility than the hadron resonance gas model, but requires input information on the quasiparticles. 
The channel treated here illustrates the importance of having guidance from an effective theory in 
reconstructing the gross features of the spectral function. In doing so, applying the Backus-Gilbert 


^ The current is here normalized as '/2V^(x), see Eq. (|^. 
® Here /s(p) = — 1)“^ is the Bose distribution. 
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method in a first step can be useful in narrowing down the region of frequency where a specihc 
ansatz for the spectral function must be made. 
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Appendix A: Tensor structure of the axial current two-point functions 


We work in the Euclidean field theory and define the correlation function in momentum space 


as 




(Al) 


The unit vector e points in the direction that defines the thermal boundary condition. It is 
e = (1,0,0,0) in the rest frame of the thermal system. By doing the change of integration variables 
X —)• —X and using translation invariance 


and symmetry under the 0(4) rotation x —)■ —x, 

{A^{-x)Au{0))e = {A^{x)Ai,{0))_i, 

we have the symmetries 

uUe, k) = (e, -k) = (-6^ -k). 


(A2) 

(A3) 

(A4) 


To write down the tensor decomposition, we have the building blocks k^ and at our disposal. 

{^fikp k^Cu'), (A5) 


We can write down four structures that respect the symmetries (A4), 

k^ku e ■ k 

We can form one projector to the subspace orthogonal to both e and k, 


/^T,t _ e_ \ _ 

'-'liv —Ofiu (p ^ i.\2 


1 - (e • fe)2/A:2 

and one projector onto the component of e which is orthogonal to k, 


“ "p , (A6) 


gT.! = a _ ^ . _ 

k^ 1 - (e.fc)2/P^^^ 


(e • k)kf, 
k^ 


eu - 


(e ■ k)k^ 
k^ 


(A7) 
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Two possible non-transverse combinations are 


„L,1 _ 






(1 \2 kfiki/ \ 

(^e^e. - (e •/c) 


(A8) 


The first one is the projector onto the direction of The second tensor, while not a projector, 
has the properties 


_ 1 


^T,t _ ^L,l _ Q 

^fia ^aiy ^fia ^au 


In summary, we can write 


k) = cj’* nT.t + cj’i nT.i + n^.i + cf, 

The argument of the C’s is (e, k), while the argument of the form factors 11 is (e • k, k"^). 


(A9) 

(AlO) 


It is helpful to be able to invert the relation Eq. (AlO) in order to project out the form factors 
individually. We find 


ttL,i _ kfj,ki, ^ 

/j,2 


= 


nT.i = 


1 




e • k ^ 

1 

1 - (e • kY/k"^ 

21 ^ 


e-k 

k2 

- (1 + (e • k f/k^) 


n'^>t = - 


nT,i ^ jjivi ^ jjL 


m 

}■ 


(All) 

(A12) 

(A13) 

(A14) 


1. Special kinematics 


When (e • k)'^ = kP', corresponding to vanishing spatial momentum in the rest frame of the 
thermal system, the projectors Cjii* and Cjii* as well as become singular. Therefore we will 
define the value of the form factors in this limit by continuity. When e and k are collinear, there 
are only two independent tensor structures. 


n^r>(£, k) = (v - 




)nT(fc2) + 


A:2 


Applying the relevant projectors as in Eqs. (A11-A14), one finds that 

nT,t ^ nT,i ^ 


n^’' = El, 


= 0 . 


(A15) 


(A16) 


When e • A: = 0, corresponding to the static correlators, Cjii* becomes equal to C^. Therefore, 
in that situation the Euclidean correlator is only sensitive to the sum of the two corresponding 
form factors, (II'^ + 11^^). Eq. (A12) nonetheless provides an unambiguous definition of II'^ if. 


expressed in the rest frame, limfcg_>o is known. The latter limit, however, requires an analytic 

continuation of the Euclidean correlator. 


2. The zero-temperature limit 

At zero temperature, it is natural to parametrize the correlation function as 

= {s^ - (ai7) 
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Applying the same projectors as in Eqs. (All'A14) onto the correlation function (A17), and 
requiring that the same result be obtained in the zero-temperature limit, we obtain 


nL’'(e. A:,fc2) - 

-^nL(fc2), 

(A18) 

U^{e-k,k‘^) - 

-^0, 

(A19) 

^^’'(^ A:,fc2) - 


(A20) 

U^’\e-k,k‘^) - 

-A n'^(A;2). 

(A21) 


3. Relation to the correlators of the pseudoscalar density 


We define 


,5“%(e^ A:) = J <fx (A“(x) P'(0)), 

(A22) 


(A23) 

We note the symmetry relations 


A^{e, k) = -A^{-e, -k), 

(A24) 

[ d\ (P(x) A,(0)), = A,{e,-k), 

(A25) 


respectively from 0(4) invariance and from translation invariance. 

In|8], Eq. (A7) and (A8), taking into account Eq. ( |A25| ) and Eq. ( |A4[ ), the Ward identities 

tA (-^ 26 ) 

(A27) 


2mAfj,{i,k) = ika n^o(e, fc), 

4m^ V{e, k) = k)ka + m(VSV') 


were 


derived. Inserting our tensor decomposition of k), we find 

2mA^{e, k) = ik^ n^’^(e • k, k‘^) + i{e ■ k) ^^^ 

4m^ V{e, k) = k"^ n'^’*(e • k, /c^) -|- m (V'V’)- 


(A28) 

(A29) 


Appendix B: On the residue of the pion pole 

In this appendix, we use the general results of the previous section in the rest frame of the ther¬ 
mal system, = (1,0,0,0). The form factors are thus functions of /cq and k"^ and the dependence 
on e is no longer indicated explicitly. All expressions for correlation functions in this section refer 
exclusively to the pion contribution. 

In ID, it was shown that the residue of the pion pole in the two-point function of Aq at vanishing 
spatial momentum is Res(wo) = In order to determine the form of the residue at finite 

momentum, we parametrize the residue as 

Res(u;k) = + Ak^). (Bl) 

To determine the parameter A, we will exploit the fact that the spectral representation of the two- 
point function of Aq in terms of real-time excitations must agree with the spectral representation 
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in terms of screening states. From the former point of view, the pion contribution to the correlator 
in momentum space takes the form 


Kik) = 


/,^(m^ + Ak^) 
kl+ujl 


(B2) 


with Wk given in Eq. (11). From the ‘screening’ point of view, the residue must be proportional to 
/cq at small kQ (here we invoke the analytic continuation in the frequency, away from the Matsubara 
values kQ = 2TTTn). This is so because the screening pion is odd under the Euclidean time reversal 
xq — t —xq, while ^0 is even. Thus we can write 

-\F\X 


^Uk) = 


+ ^ + 


(B3) 


for some parameter F to be determined. Equations (B2) and (B3) must agree when the numerators 
are evaluated at the pole, /cq = From here we learn the following. 


ii"i = 4 ’ 


A = 1. 


(B4) 


This shows in particular that the residue has the form given in Eq. (20). Essentially the same 
argument was already used in [8] to determine the residue of the pion pole in the two-point function 
of the pseudoscalar density. 




vik) = - 


1 


ko+^l 


We note that for a one-pole contribution, factorization relations such as 

\AQ{k)\^ = \v{k)\ in^oWI 


(B5) 


(B6) 


hold. The phase of can then be determined through its form at vanishing spatial momentum 
given in [8]. 


1. The pion contribution to 11^^ 


Having found the residue of the pion pole in the various two-point functions of the axial current, 
we give for completeness the pion contribution to the form factors defined in Eqs. (A11-A14), 


n^\ko,k‘^) 


Il^{ko,k‘^) 


n^\ko,k‘^) 

n^'\ko,k^) 


k^{kl+u:iy 

- V?) 

k‘^{kl+ujl) 

kl+^l 


0 . 


(B7) 

(B8) 

(B9) 

(BIO) 


The first is obtained from Eq. ( A29[), th en the second from Eq. (A28), the third by using Eq. (A13) 
and the first two results. Via Eq. (AlO), the form factors allow one to obtain the entire tensor H^. 


These calculations could be greatly expedited by using an effective Lagrangian, as written down 
in HU. However it is also instructive to derive the results above directly within QCD. 
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Appendix C: Lattice correlation functions 


xo/a 

n = 0 

n = 1 

n = 2 

n = 3 

n = 4 

n = 5 

0 

2.9921(9) X 10^ 

2.9962(9) X 10° 

3.0079(9) X 10° 

3.0260(9) X 10° 

3.0496(9) X 10° 

3.0774(9) X 10° 

1 

9.200(5) X 10^ 

9.335(5) X lOi 

9.698(5) X lOi 

1.0213(5) X 10^ 

1.0813(5) X 102 

1.1447(5) X 102 

2 

5.65(1) X 10° 

6.27(1) X 10° 

7.83(1) X 10° 

9.77(1) X 10° 

1.171(1) X 10° 

1.343(1) X lOi 

3 

1.566(6) X 10° 

1.903(5) X 10° 

2.681(5) X 10° 

3.527(4) X 10° 

4.222(4) X 10° 

4.699(5) X 10° 

4 

5.19(5) X 10-1 

7.10(4) X 10-1 

1.117(3) X 10° 

1.491(3) X 10° 

1.721(3) X 10° 

1.801(3) X 10° 

5 

4.02(5) X 10-1 

5.08(4) X 10-1 

7.11(3) X 10-1 

8.58(2) X 10-1 

9.00(2) X 10-1 

8.59(2) X 10-1 

6 

3.84(5) X 10-1 

4.34(4) X 10-1 

5.16(3) X 10-1 

5.46(2) X 10-1 

5.14(1) X 10-1 

4.44(1) X 10-1 

7 

3.77(5) X 10-1 

3.90(4) X 10-1 

3.95(3) X 10-1 

3.64(2) X 10-1 

3.06(1) X 10-1 

2.386(7) X 10-1 

8 

3.72(5) X 10-1 

3.59(4) X 10-1 

3.16(2) X 10-1 

2.54(1) X 10-1 

1.887(9) X 10-1 

1.319(6) X 10-1 

9 

3.68(5) X 10-1 

3.37(4) X 10-1 

2.64(2) X 10-1 

1.85(1) X 10-1 

1.217(8) X 10-1 

7.55(5) X 10-2 

10 

3.66(5) X 10-1 

3.23(4) X 10-1 

2.30(2) X 10-1 

1.44(1) X 10-1 

8.35(8) X 10-2 

4.59(5) X 10-2 

11 

3.65(5) X 10-1 

3.15(4) X 10-1 

2.11(2) X 10-1 

1.21(1) X 10-1 

6.37(7) X 10-2 

3.11(5) X 10-2 

12 

3.65(5) X 10-1 

3.13(4) X 10-1 

2.06(2) X 10-1 

1.14(1) X 10-1 

5.77(8) X 10-2 

2.67(6) X 10-2 


TABLE V. Renormalized correlation function G'^(xo, T, p„)/T^. All errors quoted are statistical. 
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